Dynamical control of qubit coherence: Random versus deterministic schemes 
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We revisit the problem of switching oft unwanted phase evolution and decoherence in a single 
two-state quantum system in the light of recent results on random dynamical decoupling methods 
[L. Viola and E. Knill, Phys. Rev. Lett. 94, 060502 (2005)]. A systematic comparison with 
standard cyclic decoupling is effected for a variety of dynamical regimes, including the case of both 
semiclassical and fully quantum decoherence models. In particular, exact analytical expressions 
are derived for randomized control of decoherence from a bosonic environment. We investigate 
quantitatively control protocols based on purely deterministic, purely random, as well as hybrid 
design, and identify their relative merits and weaknesses at improving system performance. We find 
that for time-independent systems, hybrid protocols tend to perform better than pure random and 
may improve over standard asymmetric schemes, whereas random protocols can be considerably 
more stable against fiuctuations in the system parameters. Beside shedding light on the physical 
requirements underlying randomized control, our analysis further demonstrates the potential for 
explicit control settings where the latter may significantly improve over conventional schemes. 

PACS numbers; 03.67.Pp, 03.65.Yz, 05.40.Ca, 89.70.-|-c 



I. INTRODUCTION 

The design and characterization of strategies for con- 
trolling quantum dynamics is vital to a broad spectrum 
of applications within contemporary physics and engi- 
neering. These range from traditional coherent-control 
settings like high-resolution nuclear jlj, i2j and molecular 
spectroscopy to a variety of tasks motivated by the 
rapidly growing field of quantum information science . 
In particular, the ability to counteract decoherence ef- 
fects that unavoidably arise in the dynamics of a real- 
world quantum system coupled to its surrounding envi- 
ronment is a prerequisite for scalable realizations of quan- 
tum information processing (QIP), as actively pursued 
through a variety of proposed device technologies . 

Active decoupling techniques offer a conceptually sim- 
ple yet powerful control-theoretic setting for quantum- 
dynamical engineering of both closed-system (unitary) 
and open-system (non-unitary) evolutions. Inspired by 
the idea of coherent averaging of interactions by means 
of tailored pulse sequences in nuclear magnetic resonance 
(NMR) spectroscopy |^, decouphng protocols consist 
of repetitive sequences of control operations (typically 
drawn from a finite repertoire), whose net effect is to co- 
herently modify the natural target dynamics to a desired 
one. In practice, a critical decoupling task is the selec- 
tive removal of unwanted couplings between subsystems 
of a fully or partially controllable composite quantum 
system. Historically, a prototype example is the elimi- 
nation of unwanted phase evolution in interacting spin 
systems via trains of 7r-pulses (the so-called Hahn-echo 
and Carr-Purcell sequences 0,11). For open quantum 
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systems, this line of reasoning motivates the question of 
whether removing the coupling between the system of in- 
terest and its environment may be feasible by a control 
action restricted to the former only. Such a question was 
addressed in Pi for the paradigmatic case of a single qubit 
coupled to a bosonic reservoir, establishing the possibil- 
ity of decoherence suppression in the limit of rapid spin 
flipping via the echo sequence mentioned above. 

The study of dynamical decoupling as a general strat- 
egy for quantum coherent and error control has since 
then attracted a growing interest from the point of view 
of both model-independent decoupling design and opti- 
mization, and the application to specific physical sys- 
tems. Representative contributions include the extension 
to arbitrary finite-dimensional systems via dynamical- 
algebraic geometric and linear-algebraic 
formulations; the construction of fault-tolerant Eule- 
rian 0| and concatenated decoupling protocols fisl. a s 
well as efficient combinatorial schemes UM [j7l llSl Il9j : 
the connection with quantum Zeno physics ; proposed 
applications to the compensation of specific decoherence 
mechanisms (notably, magnetic state decoherence |2ll | 
and 1// noise ,22i, ,2^ 2Ai .26]) and/or the removal 
of unwanted evolution within trapped-ion 1271 1231 and 
solid-state quantum computing architectures [23 • These 
theoretical advances have been paralleled by steady ex- 
perimental progress. Beginning with a proof-of-principle 
demonstration of decoherence suppression in a single- 
photon polarization interferometer |3C| . dynamical de- 
coupling techniques have been implemented alone and 
in conjunction with quantum error correction within 
liquid-state NMR QIP ISlllsal , and have inspired charge- 
based and flux-based |34| echo experiments in super- 
conducting qubits. Recently, dynamic decoherence con- 
trol of a solid-state nuclear quadrupole qubit has been 
reported [ssl ]. 

All the formulations of dynamical decoupling men- 
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tioned so far share the feature of involving purely deter- 
ministic control actions. In the simplest setting, these are 
arbitrarily strong, effectively instantaneous rotations (so- 
called bang-bang controls) chosen from a discrete group 
G. Decoupling according to Q is then accomplished by se- 
quentially cycling the control propagator through all the 
elements oiQ. If At denotes the separation between con- 
secutive control operations, this translates into a minimal 
averaging time scale Tc = \G\At, of length proportional 
to the size \G\ of G- 

The exploration of decoupling schemes incorporating 
stochastic control actions was only recently undertaken. 
A general control-theoretic framework was introduced 
by Viola and Knill in [s^ (see also js^l), based on 
the idea of seeking faster convergence (with respect to 
an appropriately defined metric) by randomly sampling 
rather than systematically implementing control opera- 
tions from G- Based on general lower bounds for pure- 
state error probabilities, the analysis of indicated 
that random schemes could outperform their cyclic coun- 
terpart in situations where a large number of elementary 
control operations is required or, even for small control 
groups, when the interactions to be removed vary them- 
selves in time over time scales long compared to At but 
short compared to T^. Furthermore, it also suggested 
that advantageous features of pure cyclic and random 
methods could be enhanced by appropriately merging 
protocols within a hybrid design. The usefulness of ran- 
domization in the context of actively suppressing coher- 
ent errors due to residual static interactions was mean- 
while independently demonstrated by the so-called Pauli 
Random Error Correction Method (PAREC), followed by 
the more recent Embedded Dynamical Decoupling Method 
- both due to Kern and coworkers • Both proto- 

cols may be conceptually understood as following from 
randomization over the Pauli group Gp = o'x, o'y, az}, 
used alone or, respectively, in conjunction with a second 
set of deterministic control operations. 

Our goal in this work is twofold: first, to develop a 
quantitative understanding of typical randomized con- 
trol performance for both coherent and decoherent phase 
errors, beginning from the simplest scenario of a single 
qubit already investigated in detail in the deterministic 
case ,^J; second, to clarify the physical picture underlying 
random control, by devoting, in particular, special atten- 
tion to elucidate the control action and requirements in 
rotating frames associated with different dynamical rep- 
resentations. The fact that the controlled dynamics re- 
mains exactly solvable in the bang-bang (BB) limit makes 
the single-qubit pure-dephasing setting an ideal test-bed 
for these purposes. From a general standpoint, since spin- 
flip decoupling corresponds to averaging over the smallest 
(nontrivial) group 22 = {0,1}, with Tc = 2At [SIlll, 
this system is not yet expected to show the full advantage 
of the random approach. Remarkably, however, control 
scenarios can still be identified, where randomized pro- 
tocols indeed represent the most suitable choice. 

The content of the paper is organized as follows. Af- 



ter laying out the relevant system and control settings 
in Sect. II, we begin the comparison between cyclic and 
randomized protocols by studying the task of phase refo- 
cusing in a qubit evolving unitarily in Sect. III. Control 
of decoherence from purely dephasing semiclassical and 
quantum environments is investigated in the main part 
of the paper. Sects. IV and V. We focus on the relevant 
situations of decoherence due to random telegraph noise 
and to a fully quantum bosonic bath, respectively. Both 
exact analytical and numerical results for the controlled 
decoherence process are presented in the latter case. We 
summarize our results and discuss their significance from 
the broader perspective of constructively exploiting ran- 
domness in physical systems in Sect. VI, by also pointing 
to some directions for future research. Additional techni- 
cal considerations are included in a separate Appendix. 

II. SINGLE-QUBIT QUANTUM-CONTROL 
SETTINGS 

Our target system 5 is a single qubit, living on a state 
space Hs — C^. The influence of the surrounding en- 
vironment may be formally accounted for by two main 
modifications to the isolated qubit dynamics. First, S 
may couple to effectively classical degrees of freedom, 
whose net effect may be modeled through a determinis- 
tic or random time-dependent modification of the system 
parameters. Additionally, S may couple to a quantum en- 
vironment E, that is, a second quantum system defined 
on a state space He with which S may become entan- 
gled in the course of the evolution. For the present pur- 
poses, E will be schematized as a bosonic reservoir con- 
sisting of independent harmonic modes. Let 1s,e denote 
the identity operator on Hs,e, respectively. Throughout 
the paper, we will consider different dynamical scenarios, 
corresponding to special cases of the following total drift 
Hamiltonian on Hs 'E)'He- 

Hoit) = Hs{t)(g,lE + ts(g>HE + HsEit), (1) 

where 

' Hs{t) = , 

< HE = Y.k^kbibk, (2) 

^ HsE{t)^ t^Oz®Y.k [9k{t)bl+gl{t)bk) . 

Here, we set ft = 1, and {i = x, y, z), 6^ and bk denote 
Pauli spin matrices, and canonical creation and annihi- 
lation bosonic operators of the kth environmental mode 
with frequency ujk, respectively, ujoit) and gk{t) are real 
and complex functions that account for an effectively 
time-dependent frequency of the system and its coupling 
to the kth reservoir mode, respectively. We shall write 

LOo{t) = UJo+ SujQ{t) , 

gkit) = gk + 6gk{t) , (3) 
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for an appropriate choice of central values wg, gk and 
modulation functions Sloq, 6gk, respectively. The adi- 
mensional parameter fi is introduced for notational con- 
venience, allowing to include (/i — 1) or not (/i — 0) 
the coupling to E as desired. Physically, because Hs{t) 
and Hssit) commute at all times, the above Hamilto- 
nian describes a purely decohering coupling between S 
and E, which does not entail energy exchange. While in 
general dissipation might also occur, focusing on pure de- 
coherence is typically justified for sufhciently short time 
scales 113, 01 and, as we shall see, has the advantage of 
making exact solutions available as benchmarks. 

Control is introduced by adjoining a classical controller 
acting on S, that is by adding a time-dependent term to 
the above target Hamiltonian, 

Ho{t) ^ Ho{t) + H,{t) ® 1e . (4) 

In our case, Hc{t) will be designed so as to implement 
appropriate sequences of BB pulses. This may be accom- 
plished by starting from a rotating radiofrequency field 
(or, upon invoking the rotating wave approximation, by a 
linearly-polarized oscillating field), described by the fol- 
lowing amplitude- and phase-modulated Hamiltonian: 



with 

vi=){t) = vm-t,)-Q{t-t,-T)]. 

Here, 0(-) denotes the Heaviside step function (defined 
as &{x) = for a; < and 6 = 1 for a; > 0), ^ and r 
are positive parameters, and tj denotes the instants at 
which the pulses are applied. If the carrier frequency is 
tuned on resonance with the central frequency, uj = loq^ 
and the phase Lpj (t) = —tootj for each j, the above Hamil- 
tonian schematizes a train of identical control pulses of 
amplitude V and duration t in the physical frame. Un- 
der the BB requirement of impulsive switching (r — > 0) 
with unbounded strength {V oo), it is legitimate 
to neglect 77o(i) (including possible ofT-resonant effects) 
within each pulse, effectively leading to qubit rotations 
about the ai-axis. In particular, a tt rotation corresponds 
to 2Vt = ±n (see also Appendix A). 

In what follows, we shall focus on using trains of BB 
TT-pulses to effectively achieve a net evolution character- 
ized by the identity operator (the so-called no-op gate). 
This requires averaging unwanted (coherent or decoher- 
ent) (Tz evolution generated by either Hs{t) or Hssit) 
or both, by subjecting the system to repeated spin- 
flips. In group-theoretic terms such protocols have, as 
mentioned, a transparent interpretation as implement- 
ing an average over the group Z2, represented on Tis as 
Q = {ge} = {IjCa;} |lOj |. The quantum operation effect- 
ing such group averaging is the projector Hg on the space 
of operators commuting with Q, leading to 

1^1 «ee ^ 



Essentially, in cyclic decoupling schemes based on Q the 
above symmetrization is accomplished through a time av- 
erage of the effective Hamiltonian determining the evo- 
lution over a cycle, Tc; in random schemes, it emerges 
from an ensemble average over different control histories, 
taken with respect to the uniform probability measure 
over Q [3^ . Neither deterministic nor stochastic se- 
quences of TT pulses achieve an exact implementation of 
Hg for a fully generic Hamiltonian as in Eqs. Ql-Q, ex- 
cept in the ideal limit of arbitrarily fast control where the 
separation between pulses approaches zero. Therefore, it 
makes sense to compare the performance attainable by 
different control sequences for realistic control rates. In 
this paper we shall focus on the following options. 

• Asymmetric cyclic protocol (A), Fig.^a). This 
is the protocol used in j^], corresponding to repeated 
spin-echoes. Cyclicity is ensured by subjecting the sys- 
tem to an even number of equally spaced 7r-pulses, ap- 
plied at tj = ^0 + J^^i i = Ij 2, . . ., in the limit r — > 0. 
The elementary cycle consists of two pulses: the first one, 
applied after the system evolved freely for an interval At, 
reverses the qubit original state and the second one, ap- 
plied a time At later, restores its original state. 

• Symmetric cyclic protocol (S), Fig. ^b). This 
protocols, which is directly inspired to the Carr-Purcell 
sequence of NMR, is obtained from (A) by rearranging 
the two TT-pulses within each cycle in such a way that the 
control propagator is symmetric with respect to the mid- 
dle point. The first pulse is applied at ti ~ to + and 
the next ones at tj =ti + (j — 1) At, with j > 1. Both the 
A- and S-protocols have a cycle time = 2At, and lead 
to the same averaging in the limit At 0. For finite At, 
however, the symmetry of the S-protocol guarantees the 
cancellation of lowest-order corrections 0(At), resulting 
in superior averaging |^. [25ll4^. 

• Long symmetric cyclic protocol (LS), Fig.^c). 
This is basically an S-protocol with a doubled control in- 
terval, At I— > 2At. Equivalently, note that this scheme 
corresponds to alternating a 7r-pulse with the identity af- 
ter every At. The cycle time becomes Tc = 4 At. For this 
amount of time, twice as many pulses would be used by 
protocols (A, S). Still, in certain cases, the LS-protocol 
performs better than the A-protocol (see Sect. V.D), 
which motivates its separate consideration here. 

• Naive random protocol (R), Fig. ^d). Ran- 
dom decoupling is no longer cyclic, meaning that the 
control propagator does not necessarily effect a closed 
path (see also for a discussion of acyclic deterministic 
schemes). The simplest random protocol in our setting 
corresponds to having, at each time tj = to -t- jAt, an 
equal probability of rotating or not the qubit that is, at 
every tj the control action has a 50% chance of being a 
TT-pulse and a 50% chance of being the identity. In order 
not to single out the first control slot, it is convenient to 
explicitly allow the value j = Q (equivalently, to consider 
a fictitious pulse Pq = t m the A-, S- and LS-protocols). 
For pure phase errors as considered, such a protocol may 
be interpreted as a simplified PAREC scheme [s^ . While 
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Deterministic Protocols 

1st cycle — 2nd cycle — - 

I I I I 
I I I I asymmetric (A) 

1st cycle — - . — 2nd cycle — • 

III 
I I I symmetric (S) 

' 1st cycle " 

(c) 

I I I I I long symmetric (LS) 

to t„+At t„+2At t„+3At t„+4At 

Random Protocols 

(d) 

I I I I I random (R) 

I I I I I 

I I I I I (e) 

I I I I I hybrid (H) 

to t„+At t„+2At to+3At t„+4At 

FIG. 1: Pictorial representation of relevant control protocols 
used for coherence control. Deterministic pulses are indicated 
with full vertical lines, while random pulses correspond to 
dashed vertical lines. 

we will mostly focus on this naive choice in our discussion 
here, several variants of this protocol may be interest- 
ing in principle, including unbalanced pulse probabilities 
and/or correlations between control operations. 

• Hybrid protocol (H), Fig. ^e). Interesting con- 
trol scenarios arise by combining deterministic and ran- 
dom design. The simplest option, which we call "hy- 
brid" protocol here, consists of alternating, after every 
At, a TT-pulse with a random pulse, instead of the iden- 
tity as in the LS-protocol. For our system, in the em- 
bedded decoupling language of [s^ , this may be thought 
of as nesting the A- and R-protocols. In group-theoretic 
terms, the H-protocol may be understood as randomiza- 
tion over cycles |36| . A complete asymmetric cycle may 
be constructed in two ways, say Ai and A2. Cycle Ai 
corresponds to traversing Q in the order {l,ax) that is, 
free evolution for Ai; first pulse; free evolution for At; 
second pulse - the cycle being completed right after the 
second pulse. Cycle A2 corresponds to the reverse group 
path, (ffx, l).Thus, we have: pulse; free evolution for At; 
second pulse; and another free evolution for At - the sys- 
tem should be observed at this moment before any other 



pulse. The H-protocol consists of uniformly picking at 
random one of the two cycles at every instant t2j, where 
J =0,1,... 



III. RANDOMIZED PHASE REFOCUSING IN 
AN ISOLATED QUBIT 

A single-qubit evolving according to unitary dynamics 
(/i = in Eq. |5J) provides a pedagogical yet illustrative 
setting to study dynamical control. Since the goal here is 
to refocus the underlying phase evolution, the analysis of 
this system provides a transparent picture for the differ- 
ences associated with deterministic and random pulses. 
It also simplifies the comprehension of the results for the 
more interesting case of a single-qubit interacting with a 
decohering semiclassical or quantum environment, where 
the control purpose becomes twofold: phase refocusing 
and decoherence suppression. 

A. Time-independent qubit Hamiltonian 

We begin by considering the standard case of a time- 
independent target dynamics, ujQ{t) = loq for all t. For 
all the control protocols illustrated above, the system 
evolves freely between pulses, with the propagator 

Ua{tj+i,tj) = e-"^o{t,+i-t,)a,/2 ^ 

whereas, during a pulse, it is only affected by the control 
Hamiltonian. The propagator for an instantaneous pulse 
applied at time t — tj will be indicated by Pj . Let 

Pit)^ E Pimimiml (6) 

denote the qubit density operator in the computational 
basis {|0),|1)}, with a^lO) = |0) and a^\l) = 
The relevant phase information is contained in the off- 
diagonal matrix element poi{t). If p(io) is the initial 
qubit state, the time evolution after M control intervals 
under either deterministic or randomized protocols, 

p{tM) = U{tM,to)pito)U\tM,to) , (7) 

is dictated by a propagator of the form 



U{tM,to) = T exp |-i ^ [Ho + H,iu)] rfwj = PnUoitM ,tM-i)PM-iUa{tM-utM-2) ■ . . PiUoiti,to)Po 

= {PmPm-1 ■ . ■ PiPo){Pm-i . . . PiPo)^Uo{tM,tM-i)(PM-i ■ ■ ■ PiPa) ■ ■ ■ Uo(t2,ti){PiPo)P^Ua{ti,to)Po , (8) 

V ' 

Uc{tM) 

I 

where T indicates, as usual, time ordering. Recall the basic idea of deterministic phase refocusing. 

For the A-protocol, Pq = 1 and Pj ~ exp{—iTrax/2), 
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j > 1 [see Appendix A] . Exact averaging is then ensured 
after a single control cycle, thanks to the property 

Thus, the total phase that the qubit would accumulate in 
the absence of control is fully compensated, provided that 
N complete cycles are effected (that is, an even number 
M = 2N of spin flips is applied). The overall evolution 
implements a stroboscopic no-op gate. U(tM,to) = 1, 
(tM - to) = MAt = NTc, as desired Notice that 

the identity operator is also recovered with the S- and LS- 
protocols after their corresponding cycle is completed. 

1. Control performance in the logical frame 

In preparation to the randomized protocols (R,H), it is 
instructive to look at the system dynamics in a different 
frame. In particular, a formulation which is inspired by 
NMR 0] is the so-called toggling-frame or logical-frame 
picture, which corresponds to a time-dependent interac- 
tion representation with respect to the applied control 
Hamiltonian. Let 

Uc{t,to) =T exp|-i^ Hciu)du^ (10) 

denote the control propagator associated to Hc{t). Then 
the transformed state is defined as 

p(1^ = Ul{tM)p{t)Uc{tM), (11) 

with tilde indicating henceforth logical-frame quantities. 
At the initial time the two frames coincide and 
p(io) — p{to)- The evolution operator in the logical frame 
is immediately obtained from Eqs. ([TJl and Hll|l . 

Uit,tn) ^uUt,to)U{t,to) , (12) 

with 

U{t,to) ^ T cxp l^-i [Ul{u)HoU,{u)]duj . (13) 

That is, the control field is explicitly removed from the 
effective logical Hamiltonian. Because, for BB multipulse 
control, 

Uc{tM,to) - PmPm-1 ■ ■ ■ PiPn , (14) 

the expression for the logical frame propagator may sim- 
ply be read off Eq. (jSJ , yielding 

UitM,to) - T n 'Pluit,+i,t,)V, j , (15) 

in terms of the composite rotations 

V, = P,P,-i ...PiPo, J = 0, . . . , M - 1 . 



For cyclic protocols, Uc{tM,to) = exp{—iMiTax/'2) = 
1 {M even) that is, the logical and physical frames over- 
lap stroboscopically in time. Thus, U{tM,ta) = 1 and 
phase refocusing in the logical frame is equivalent to phase 
refocusing in the physical frame. 

Now consider the evolution under the randomized pro- 
tocols. The first pulse occurs at to, so after a time interval 
tM—to has elapsed, M-\-l pulses have been applied. Since 
the final goal is to compare random with cyclic controls, 
we shall take M even henceforth. At time t = tM, pop- 
ulation inversion may have happened in general in the 
physical frame. This makes it both convenient and nat- 
ural to consider the logical frame, where inversion does 
not happen, as the primary frame for control design. The 
evolution operator in this frame may be expressed, using 
Eq. lO, as 

U{tM,to) = expl -i'^^^^a^^ Xjl , (16) 

j=o J 

where 

X, = (-1)^°+^^+-+^^ j=0,...,M-l, (17) 

is a Bernoulli random variable which accounts for the 
history of spin flips up to tj in a given realization. For 
each m = 1, . . . , j, if a spin flip occurs at time tm, then 
A„i — 1 and Pm — —i<Jx, otherwise A,„ = and P„i = 1. 
Equivalently, Xj will take the values -1-1 or —1 with equal 
probability, depending on whether the composite pulse 
Vj is the identity or a 7r-pulse. 

Let k be an index labelling different control realiza- 
tions. For a fixed k, the qubit coherence in the logical 
frame is given by 

f^itM) = e-'"°^*S"o''^^"poi(io) . (18) 

This expression provides the starting point for analyzing 
control performance. For the A-protocol, the only pos- 
sible realization has Xj — (^1)"* and leads to the trivial 
result poi(iM) — Poi{to)- For the R-protocol, realiza- 
tions corresponding to different strings of A's filling up 
M places give, in general, different phases and an en- 
semble average should be considered. If the statistical 
ensemble is large enough, the average performance may 
be approximated by the expected performance^ which is 
obtained by averaging over all possible control realiza- 
tions and will be denoted by E( ). The calculation of the 
expectation value is straightforward in the unbiased set- 
ting considered here. Since, for each realization, Xj ~ +1 
or —1 independently of the value of its predecessor Xj-i, 
the following expression is found: 

E(poi(iM)) -poi(io)[cos(woAi)]^. (19) 

Several remarks are in order. Under random pulses, 
the phase accumulated during the interval t^—tg — MAt 
is, on average, completely removed, regardless of the At 
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value. An important distinction with respect to the de- 
terministic controls, however, is that now the different 
phase factors carried by each stochastic evolution may 
interfere among themselves, causing the ensemble aver- 
age to introduce an effective phase damping. In general, 
let us write the ensemble expectation in the form 



E 



(poi(iM)) 
Poi(^o) 



(tm ,ta) —r,{tM ,to) 



(20) 



for real functions (/)*(t), r*(t) [i^. Here, 4>-t{tM,to) = 
0, whereas r»(tM,io) = — Mln[cos(a;oAt)]. Complete 
dcphasing occurs when ujoAt = Z'7r/2, with V odd, while 
for loqAI — In, with I g Z, r*(t) = 0. Whenever exact 
knowledge of the frequency loq and precise control over 
the time interval At are available, the R-protocol can 
be made to achieve exact averaging, like the A-protocol, 
under the additional synchronization condition that 

At = In/ujQ , I eZ. 

In situations where such a synchronization is not eas- 
ily accessible, one may still look for a general condition 
under which the R-protocol avoids ensemble dephasing. 
Taking a Taylor expansion of Eq. (|19|l yields 

Lj^itM - to)At = Lo^itM ~ tof/M < 1 . (21) 

In principle, this requirement may be fulfilled by making 
tui and/or At sufficiently small. Interestingly, the condi- 
tion of Eq. (|21|) is directly related to the bound obtained 
in Theorem 1 of ^] for the worst-case pure-state error 
probability, defined by 

et = max{et{m} = 1 - min E(Tr(p(to)p(i))) , (22) 

where the latter term is the usual input-output state fi- 
delity 0. In the limit where \\Hoit)\\ltAt 1, where 
I2 — max |eig(yl)|, = A\ Theorem 1 implies 



et=o(\\Ho{t)\\ltAty (23) 
On the other hand, using Eq. H19|l we obtain 



et„(|V)) = 2 poo(io)pii(io) - \poiito)\^[cos{uoAt)] 



M 



For Lij'^{tM — tQ)At <^ 1, the above expression gives 

£tM(IV^» ~ |Poi(io)|'^o(iM - to) At , (24) 
et,,=oLl{tM-to)At) , (25) 



which makes the connection with Eq. H21|l manifest. 

It remains to discuss the performance of the H- 
protocol. The freedom of not always effecting a spin flip 
after every At, which is one of the appealing features 
of the R-protocol, is still partially present here. On the 
other hand, since a spin flip does occur at every tm with 



m odd, which leads to Xj = ~Xi-i for J odd, any re- 
alization of this protocol completely refocuses the qubit 
[see Eq. CHU, so (j)*{tM,to) = and T^{tM,to) = 0. 
Accordingly, in the logical frame, the H-protocol is opti- 
mal, combining the absence of phase damping of cyclic 
schemes with the flexibility of random pulses. 



2. Ensemble averages: General remarks 

In practice, we deal with the average performance of 
a statistical ensemble of size K. To evaluate the sample 
size that guarantees a desired margin of error 5 [4^ , we 
invoke the central limit theorem. Because different real- 
izations are independent, the latter ensures that the av- 
erage performance is distributed normally with a mean 
value equal to the expected performance and standard 
deviation given by a/\/K, where a is the standard devi- 
ation for all realizations. Thus, if we want, with proba- 
bility (1 — e), that the average performance differs from 
the expected performance by no more than 6, the sample 
size must be at least as large as 



o 



(26) 



where is the value of the standard normal variable 
which has a probability e/2 of being exceeded. Taking a 
Taylor expansion of Eqs. (|18l) and p9|l . we can show that 



a = Oiujoy/{tM -ta)At) for Lu^{tM-ta)At<^l. 

Thus, the number of realizations required to ensure a 
specified degree of precision decreases as At. 

It is interesting to observe that the ensemble average 
may be interpreted as effecting a quantum operation. 



^(poiitAl)) 



p{to) 



with 



^[7(fe)t[7(fc) = i 



and random coefficients a^'^^ , P^''^ which may be derived 
from Eq. 



3. Control performances in the physical frame 

Finally, it is important to compare the average coher- 
ence element in the logical and physical frames. Dephas- 
ing is a more delicate issue in the Schrodinger picture, 
because spin population is not necessarily conserved and 
PQi{tM) rnay be related to poi(^o) or to pio{to), depend- 
ing on how many 7r-pulses occur. If, after an interval 
tM — to, an even number of spin-flips have happened, we 
recover Uc{tM,to) = 1 as in the cyclic case, but an odd 
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number of flips leads instead to UdtM^to) = ±iax- By 
recalling Eq. 112|l . for randomized schemes we find 



E(poi(^m)) = 



(27) 



where T^{tM,to) = for the H-protocol. Thus, the agree- 
ment between the expected results in the two frames de- 
pends on the initial qubit state. Results are identical if 
Poi(^o) is real, but differ otherwise. The worst scenario 
occurs if poi(^o) is purely imaginary, as the average in the 
physical frame vanishes. This reflects the fact that the 
net evolution may be represented by a quantum operation 
which flips the state of the qubit with 50% probability, 
and leaves it alone otherwise. Clearly, knowledge of the 
control history allows the system to be deterministically 
returned in the physical frame for any realization, if de- 
sired. That is, having a classical register that records 
the total number of spin flips may be used to select re- 
alizations that guarantee a good performance of random 
pulses also in the physical frame for any initial state. For 
example, if only realizations with an even number of spin- 
flips are selected, the results in both frames are equal, 
Hpoi{tM)\Uc{tM,to = 1) = E{poi{tMj), as desired. 

To summarize: In the logical frame, refocusing the un- 
wanted phase evolution is possible with any of the pro- 
tocols we considered. The R-protocol, however, intro- 
duces an average ensemble dephasing, which may only be 
prevented by precisely tuning At — In/uJo, with ^ S Z, 
or by assuring that At <C l/[wo(ij\/ — ^o)]- This im- 
plies the appearance of a time scale requirement which 
is not present when dealing with deterministic controls, 
nor with the H-protocol. In the physical frame, state- 
independent conclusions regarding the system behavior 
may be drawn conditionally to speciflc subsets of con- 
trol realizations. Overall, the H-protocol emerges as an 
alternative of intermediate performance, which partially 
combines advantages from determinism and randomness. 



B. Time-dependent qubit Hamiltonian 

We now consider the more interesting case where the 
qubit frequency is time dependent, a;o(i) = ij^o + S^oi't)^ 
SuJo{t) = u}oG{t) being a deterministic (but potentially 
unknown) function. This could result, for example, from 
uncontrolled drifts in the experimental apparatus. 

While all protocols become essentially equivalent in the 
limit AI oo, searching for the best protocol becomes 
meaningful in practical situations where pulsing rates are 
necessarily finite. Under these conditions, the determin- 
istic protocols described so far will no longer be able, 
in general, to completely refocus the qubit. This would 
require a very specific sequence of spin fiips for each par- 
ticular function duJo{t), which would be hard to construct 
under limited knowledge about the latter. On the other 
hand, the average over random realizations does remove 
the phase accumulated for any function Sujoit), making 
randomized protocols ideal choices for phase refocusing. 



As a drawback, however, ensemble dephasing may be in- 
troduced. Thus, the selection of a given protocol will be 
ultimately dictated by the resulting tradeoffs. 
The propagator in the logical frame now reads 



U(tM,to) = exp 



t — a 
2 



3=0 



tj + i 
(1 



G{u))du 



(28) 

which reduces to Eq. (|16|l when G{t) = 0. 

Some assumptions on both the amplitude and fre- 
quency behavior of G(t) are needed in order to draw some 
general qualitative conclusions. First, if \G{t)\ ^ 1, the 
analysis developed in the previous section will still ap- 
proximately hold. In the spirit of regarding luq as a cen- 
tral frequency, we will also discard the limit \G{t)\ ^ 1, 
and restrict our analysis to cases where maxf \G{t)\ ~ 1. 
If G{t) is dominated by frequency components which are 
very fast compared to tq ~ ujq^, the effect of G{t) may ef- 
fectively self-average out over a time interval of the order 
or longer than uJq^ . In the opposite limit, where the time 
dependence of G{t) is significantly slower than wq, deter- 
ministic controls are expected to be most efficient in refo- 
cusing the qubit, improving steadily as At decreases. In 
intermediate situations, however, the deterministic per- 
formance may become unexpectedly poor for certain, in 
principle, unknown values of At. These features may be 
illustrated with a simple periodic dependence. Suppose, 
for example, that G{t) — sm{pLLjot), and p g R. For a 
fixed time interval tf~tQ < 7t/{pujo), a significant reduc- 
tion of the accumulated phase is already possible with 
few deterministic pulses. However, care must be taken 
to avoid unintended "resonances" between the natural 
and the induced sign change. For the A-protocol, this ef- 
fect is worst at At = tt/{pujq), in which case the control 
pulses exactly occur at the moment the function changes 
sign itself, hence precluding any cancellation of G{t). 

With the R-protocol, ensemble dephasing becomes the 
downside to face. The ensemble average now becomes 



M-l 

n 



cos < LUq 



At^ 



G{u)du 



(29) 



In the absence of time dependence, phase damping is 
minimized as long as Eq. (|21|l holds. Under the above 
assumptions on G{t), the condition remains essentially 
unchanged, in agreement with the fact that the accuracy 
of random averaging only depends on ||i?o(0ll2 

Refocusing is also totally achieved with the H-protocol. 
However, unlike in the case of the R-protocol, the ensem- 
ble average no longer depends on the time independent 
part of the Hamiltonian, but only on the function G(t), 
making the identification of precise requirements on At 
harder in the absence of detailed information on the lat- 
ter. We have 



A/-2 

-r.(tM,to) ^ -Q cos<; Wo 

i=0,2,4,... 



G(u)du 



(30) 
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Fig. 12 illustrates the points discussed so far. The 
sinusoidal example is considered, and we contrast the 
two aspects to be examined: the top panels show the 
phase magnitude |0*(tA/, io)|i which is optimally elimi- 
nated with random pulses, while the bottom ones give 
the dephasing rate e^'"*'-**''*"\ which is inexistent for 
deterministic controls. The interval between pulses is 
fixed. At — l/(10a;o)i and the protocols are compared 
for two arbitrary, but relatively close values of the os- 
cillation frequency rate: p = 20^/2 and p = IOtt. The 
deterministic control is very sensitive to slight changes 
of the drift and at certain instants may behave worse 
than if pulses were completely avoided. Similarly, the H- 
protocol, even though more effective than the R-protocol 
in this example, also suffers from uncertainties related 
to G{t). On the contrary, deviations in the performance 
of the R-protocol are practically unnoticeable, making it 
more robust against variations in the system parameters. 




FIG. 2: (color online) Accumulated phase (upper panels) and 
dephasing rate (lower panels) in the absence of control [solid 
(green) line], and under the A- [(blue) stars], R- [(black) cir- 
cles], and H- [(purple) plus] protocols in the logical frame, for 
G{t) = sin(pwoi) and At = l/(10a;o). Left panels: p = 20^; 
right panels: p = IOtt. Average taken over 10'^ realizations. 
In this and all simulations that will follow, we set to = 0. 

As a further illustrative example, we consider in Fig.|2| 
the following time dependence for the qubit: 

uj^{t)=uj^[l + G{t)]D{t). (31) 

The left panels have, as before, D{t) — 1, while for the 
right panels 

^(t) = (-1)L10"«*/3J. (32) 

A fixed time tf — 2/uj() is now divided into an increas- 
ing number M of intervals At. Here, selecting the most 
appropriate protocol depends on our priorities concern- 
ing refocusing and preservation of coherence. We may, 
however, as the right upper panel indicates, encounter 
adversarial situations where the time dependence of the 



qubit frequency is such that not acting on the system is 
comparatively better than using the A-protocol. Clearly, 
depending on the underlying time dependence and the 
pulse separation, such poor performances are also ex- 
pected to occur with other deterministic protocols. In 
addition, notice that, consistent with its hybrid nature, 
the H-protocol may perform worse for values of At where 
the deterministic control becomes inefficient (compare 
right upper and lower panels). In similar situations, from 
the point of view of its enhanced stability, the R-protocol 
turns out to he the method of choice. 
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FIG. 3: (color online) Accumulated phase (upper panels) and 
dephasing rate (lower panels) in the absence of control [solid 
(green) line], and under the A- [(blue) stars], R- [(black) cir- 
cles], and H- [(purple) plus] protocols in the logical frame. 
The time interval considered is tf = 2/ajo; {(^oAt)~^ = M/2; 
and G{t) — sin{pLuot), where p — 10. The drift in the right 
panels includes D{t) = (_l) L^o^ot/sj Average taken over all 
possible realizations. 

To summarize: An isolated qubit with time-dependent 
parameters provides the simplest setting where advan- 
tages of randomization begin to be apparent, in terms of 
enhanced stability against parameter variations. On av- 
erage, phase is fully compensated, and ensemble dephas- 
ing may be kept very small for sufficiently fast control. 
Similar features will appear for a qubit interacting with 
a time- varying classical or quantum environment, as we 
shall see in Sects. IV. C and V.E. 



IV. RANDOMIZED CONTROL OF 
DECOHERENCE FROM A SEMICLASSICAL 
ENVIRONMENT 

Qubit coherence is limited by the unavoidable influence 
of noise sources. Within a semiclassical treatment, which 
provides an accurate description of decoherence dynam- 
ics whenever back-action effects from the system into the 
environment can be neglected, noise is modeled in terms 
of a classical stochastic process, effectively resulting in 
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randomly time-dependent systems. Typically, external 
noise sources, which in a fully quantum description are 
well modeled by a continuum of harmonic modes (see 
Sect. V), are represented by a Gaussian process. Here, 
we focus on localized noise sources, which may be intrin- 
sic to the physical device realizing the qubit - notably, 
localized traps or background charges, leading to a quan- 
tum discrete environment. In this case, non-Gaussian 
features become important, and are more accurately rep- 
resented in terms of noise resulting from a single or a 
collection of classical bistable fluctuators - leading to so- 
called random telegraph noise (RTN) or 1//- noise, re- 
spectively. Beside being widely encountered in a vari- 
ety of different physical phenomena 0, 0, such 
noise mechanisms play a dominant role in superconduct- 
ing Josephson-junction based implementations of quan- 
tum computers 

Recently, it has been shown that RTN and l//-noise 
may be significantly reduced by applying cyclic sequences 
of BB pulses mil mill 113. We now extend the anal- 
ysis to randomized control. As it turns out, random de- 
coupling is indeed viable and sometimes more stable than 
purely deterministic protocols. While a detailed analysis 
of randomized control of genuine 1// noise would be in- 
teresting on its own, we begin here with the case of a sin- 
gle fluctuator. This provides an accurate approximation 
for mesoscopic devices where noise is dominated by a few 
fluctuators spatially close to the system [23, ISSL [SJ, [S^ ■ 
Let the time-dependent Hamiltonian describing the noisy 
qubit be given by Eqs. Q-©, where ^ = and 



Sujoit) = RTN(i) 



(33) 



characterizes the stochastic process, randomly switching 
between two values ±v/2, v > 0. We shall in fact consider 
a semjrandom telegraph noise, that is, we assume that 
the fluctuator initial state is always +v/2. The switching 
rate from ±?;/2to=F''^/2is denoted by 7ip , with 7+ + 7- = 
7. We shall also assume for simplicity that 7+ = 7_, 
corresponding to a symmetrical process. The number 
of switching events n{t,0) in a given time interval t is 
Poisson-distributed as 
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P{n{t,0) = k} 
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FIG. 4: (color online). Decoherence rate |^(t)| = e~'"'*' from 
a symmetrical bistable fluctuator. Several values of g — v/-y 
are considered, resulting from changing the coupling strength 
V at fixed switching rate 7 = 1 a.u. 



where g < 1 {g > 1) corresponds to a fast (slow) fluctu- 
ator. Given the initial condition for the fluctuator Sp^, 
where £pg = ±1 stands for the fluctuator initially in state 
±u/2, Z{t,to) may be calculated as [5ll |. 

Z(t,io) = Ce-*(i-")(*-*'') + (l-C)e-*(i+")(*-*°\ (35) 
where 

^1 - 52 + 2ig£p C={l + a~ig£p^)/{2a), (36) 

and £p = (7- — 7+)/7 is the equilibrium population dif- 
ference. Note that for a symmetrical telegraph process, 
the only difference between the results for a fluctuator ini- 
tially in state +v/2 or —v/2 is a sign in the above phase 6. 
The decoherence rate for a slow fluctuator is much more 
significant than for a fast fluctuator. This has been dis- 
cussed in detail elsewhere ^51], and has been reproduced 
for later comparison with the controlled case in Fig. 01 
where several values of g are considered. A fast fluctua- 
tor behaves equivalently to an appropriate environment 
of harmonic oscillators and noise effects are smaller for 
smaller values of g, whereas for a slow fluctuator the de- 
coherence function saturates and becomes ^7(^ — ^0)- 



Semiclassically, dephasing results from the ensemble 
average over different noise realizations. This leads to 
the decay of the average of the coherence element. 



(P0l(^)) 

Poi(to) 



(34) 



Here, the average over RTN realizations is represented 
by ( ) and should be distinguished from the average over 
control realizations, which, as before, is denoted by E. 
The dephasing factor r(i,to) and the phase 6{t,tQ) have 
distinctive properties depending on the ratio g = u/7. 



A. Deterministic and randomized controls in the 
interaction picture 

Here we compare the reduction of RTN under the ac- 
tion of the A-, H-, and R-protocols. In order to isolate 
the effects of the noise, it is convenient to first carry out 
the analysis in the interaction picture which removes the 
free dynamics uJoaz/2. The density operator becomes 



p'{t)^U\t,to)p{to)U'\t,to) 



(37) 



where U^{t,to) = exp{iujQtaz/2) and the superscript / 
will refer to the interaction picture henceforth. The free 
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propagator between pulses is now 



U\tj+utj)^T exp. 



H^{u)du 



(38) 



with i?Q(i) = (5u;o(i)crz/2; while at tj, we have [see Ap- 
pendix A] 



exp 



exp 



'lAi —(To, 



exp 



.ujotj 



(39) 

A second canonical transformation into the logical 
frame is also considered, so that (as before) realizations 
with an even or an odd number of total spin flips are 
treated on an equal footing. We will refer to the combi- 
nation of the two transformations as the logical-IP frame. 
Similarly to Eq. ()12|l . the interaction and the logical-IP 
frame propagators are related as 

U\t,to)^Ul\t,to)U\t,to). (40) 

This leads to the following propagators at Im, 



Ul{tM, to) 
U\tM, to) 



,M-1 . 

T[X{vYu\t,^,^t,)V]^ 



3=0 



where 



V' 



pi pi 

^3 ^3-1 ■ ■ 



Our goal is to compute the ratio 
EiipUtM))) ^ 



J =0, 



,M- 1 



F{tM,to) — 



Poi(io) 



.to). 



^ {EjpiAtM))) 
Poiito) 



(41) 



where k labels, as before, different control realizations. 
Note that interchanging the order of the averages does 
not modify the results if all pulse realizations are con- 
sidered and the number of RTN realizations is large 
enough. With 10^ switch realizations no significant vari- 
ations were found by interchanging the averages. 

The decoherence rate \F{tM,to)\ for the three selected 
protocols is shown in Fig. |S1 where a time tf = IO/7 
was fixed and divided into an increasing number M of 
intervals At. The left panels are obtained for three 
slow fluctuators, g — 5,3,2, and the right panels for 
g = 1.1, 0.8, 0.1. These are the six different noise regimes 
considered in Ref. ^] , where the A-protocol was studied. 
The authors concluded that once At ^ I/7, V(tM^to) 
scales with 17^, while for At > I/7, BB pulses are still ca- 
pable of partially reducing noise due to a fast fluctuator, 
but are mostly inefficient against slow fluctuators. Here, 
we verified that among all possible realizations of pulses 
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FIG. 5: (color online) Decoherence rate from a symmetri- 
cal bistable fluctuator with g = 5,3,2,1.1,0.8,0.1. (Green) 
Solid lines represent the analytical results from IsH in the 
absence of control. (Blue) stars: A-protocol; (black) cir- 
cles: R-protocol; and (purple) plus: H-protocol. Averages 
are taken over 10^ RTN realizations and all possible 2*^ pulse 
realizations. The time interval considered is tj = IO/7, thus 
(7At)-i = A//10. 



separated by the same interval At, the realization cor- 
responding to the A-protocol yields the largest value of 
|F(tA/,io)|, whereas absence of pulses gives, as expected, 
the smallest value. This justifies why, in terms of average 
performance for finite At, we have in decreasing order: 
A-, H-, and R-protocol; while for M 00, different pro- 
tocols are expected to become equivalent. 

In terms of refocusing the unwanted phase evolu- 
tion, the above randomized protocols are optimal, since 
Arg((E(/5^i(tAf)))//3oi(io)) = 0, while the phase magni- 
tude (5(tM,io) for the A-protocol is eventually compen- 
sated as M increases. This is shown in Fig.El Notice also 
that the absolute phase is very small for fast fluctuators. 

Instead of fixing a time ty, an alternative picture of the 
performance of different protocols may also be obtained 
by fixing the number of intervals M, as in Figs. [3 (left 
panels). As expected, larger M leads to coherence preser- 
vation for longer times. Still another option is to fix the 
interval between pulses At, as in Fig. (right panel). 
As before, the A-protocol shows the best performance, 
followed by the H- and R-protocols. 



B. Randomized control in tlie physical frame 

If the interaction picture is not taken into account, 
complete refocusing is again guaranteed, on average, 
when either the R- or H-protocols are used. However, 
for the R-protocol, the qubit frequency plays a delicate 
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FIG. 6: (color online) Phase offset from a symmetrical 
bistable fluctuator with <; = 5, 3, 2, 1.1, 0.8, 0.1. (Green) Solid 
line: Analytical results in the absence of control pulses; (blue) 
stars: A-protocol; both R- and H-protocols have phase equal 
to zero. The time interval considered is tf — IO/7, so 
{'yAt)~^ — M/10. Averages computed as in Fig.|Sl 




FIG. 7: (color online) Decoherence rate from a symmetrical 
bistable fluctuator with g = 1.1. Left panels: M = 10 (top), 
M = 30 (bottom). Right panel: At = I/7. (Green) Solid 
line: Analytical results in the absence of control. (Blue) stars: 
A-protocoI; (black) circles: R-protocol; and (purple) plus: H- 
protocol. Averages are taken over 10* RTN realizations and 
10^ pulse realizations. 



role in the resulting dephasing process. We now have 

H{pQi{tM))) = Poi{to) ■ (42) 
^^-lujoAtY: fiT, ' xf +^'5<'=' (tM , to) g-rt"' (tM ,to) ) ^ 

which may be further simpUfied as follows. Among the 
2^ pulse realizations existing in the logical frame, there 
are pairs, say corresponding to labels k and k', where 

X'^o^ = 1 and A^,'''^ = 0, while xf^ = xf'^ with l<j< 
M — 1, which leads to J2j=i^ xf"* — — J2j=i^ xf ^ ■ Be- 



sides, since we are considering a semirandom telegraph 
noise, a pulse at tg is equivalent to switching the fluctua- 
tor from the initial state +v/2 to —v/2, whose net effect 
is simply a change in the sign of the phase 5(tj\/,to). 
Therefore, we may write 



E((poi(iAf)» 



Poi{to 



2M-1 

{cos[4f-i^oAi-(5W(tM,to) 



k=l 

where 



(43) 
^r<'=)(tM.to)| 



M-1 

E 



X 



(fc,l) 



M > 2, 



(44) 



andxr"^-(-l)'^''+'^''+-+'^'' 



(fc,i) ^ 
3 

In the physical frame, we find correspondingly 
Poi(io) + Pio(io) 



mpoiitm)) - 



El 

k=\ 



cos 



2 • 2*^-1 
4f^c^oA<-<5W(<M,io) 



(45) 
-r(*'(tA/,«o)| _ 



Contrary to the result obtained in the absence of noise, 
Eq. p9|l . the additional realization-dependent phase shift 
^^^^ (tMyto) now remains. While, on average, this phase is 
removed in the limit where At — *■ 0, for finite control rates 
S'''^\tM,to) may destructively interfere with the phase 
gained from the free evolution, potentially increasing the 
coherence loss. Identifying specific values of At where 
such harmful interferences may happen for the given RTN 
process is not possible, which makes the results for the 
R-protocol with finite At unpredictable in this case. 

While the above feature is a clear disadvantage, it 
is avoided by the H-protocol. For each realization, the 
phase accumulated with the free evolution is completely 
canceled, so the result in the logical-IP frame is equal to 
that in the logical frame: E((poi(iM))) = £{{poi{tM)))- 
If access to a classical register that records the total num- 
ber of spin flips is also available, this equivalence between 
frames may be further extended to the physical frame. 
Additionally, as already found in Sect. IIIB, randomized 
protocols tend to offer superior stability. 



C. Deterministic bursts of switches 

Let us illustrate the above statement through an ex- 
ample where the noisy dynamics of the system is slightly 
perturbed. Suppose that, moving back to the interaction 
picture, the noise process is now 



where 



, , , RTN(t) 



D{t) = (-l)L47t/25jL47t/5j 



(46) 



(47) 
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Physically, D{t) describes a sequence of six instanta- 
neous switches, equally separated by the interval 5/(47), 
restarting again at every instant 25k/ (ij), k being an 
odd number. This process may be viewed as bursts 
of switches of duration 25/(47) followed by an interval 
25/(47) of dormancy. The resulting behavior for g — 1.1 
in the logical- IP frame is depicted in Fig. |H1 
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FIG. 8: (color online) Decoherence rate (upper panels) and 
phase offset (lower panels )in the logical-IP frame for a single 
fiuctuator with g = 1.1 subjected to a disturbance as given in 
Eq. (|I7J. Left panels: Fixed tf = IO/7, so (7At)"^ = M/10. 
Right panels: fixed = 5/(47). (Green) Solid line: Results 
in the absence of control pulses; (blue) stars: A-protocol; 
(black) circles: R-protocol; and (purple) plus: H-protocol. 
The average phase for both R- and H-protocols is zero. Left 
panels: Averages are taken over 10^ RTN realizations and all 
possible 2*^^ pulse realizations. Right panels: Averages are 
taken over 10* RTN realizations and 10'^ pulse realizations. 



the latter case. In the logical and physical frames, 
the R-protocol, besides showing the poorest performance 
among the three considered schemes, may also lead to 
dangerous interferences between the qubit frequency and 
the phase gained from the free evolution. Such problem, 
however, does not exist for the H-protocol. The benefits 
of randomization are most clear when limited knowledge 
about the system dynamics is available and determinis- 
tic control sequences may be inefficient in avoiding un- 
wanted "resonances". Combining protocols, where we 
gain stability from randomness, but also avoid the free 
phase evolution, is desirable especially when working in 
the physical frame. In this sense, the H-protocol emerges 
as a promising compromise. 



V. RANDOMIZED CONTROL OF 
DECOHERENCE FROM A QUANTUM BOSONIC 
ENVIRONMENT 

We now analyze the case of a genuine quantum reser- 
voir, where decoherence arises from the entanglement 
between the qubit and the environment. The relevant 
Hamiltonian is given by Eqs. CQl-lO with wo(i) = "^o 
and fi = 1. In the semiclassical limit, the effects of the 
interaction with the bosonic degrees of freedom may be 
interpreted in terms of an external noise source whose 
fluctuations correspond to a Gaussian random process. 

A detailed analysis of deterministic decoherence sup- 
pression for this model was carried out in (see also [5q| 
for an early treatment of the driven spin-boson model in a 
nonresonant monochromatic field and (57l | for related dis- 
cussions of dynamically modified relaxation rates). Here, 
we discuss how randomized decoupling performs. 



With deterministic control, the rate of noise suppres- 
sion quickly improves as the separation between pulses 
shrinks (left upper panel), until a certain value. At = 
tf/8 — 5/(47), where it suddenly shows a significant re- 
coil, becoming almost as bad as simply not acting on 
the system at all. Equivalently, by fixing At — 5/(47), 
the performance of the A-protocol becomes very poor for 
t > 10/7 (right upper panel). In practice, detailed knowl- 
edge of the system dynamics might be unavailable, mak- 
ing it impossible to predict which values of At might be 
adverse. Randomized schemes, on the other hand, are by 
their own nature more stable against such interferences. 
As seen from the figure, the R-protocol shows a slower, 
but also more consistent improvement as At decreases, 
and might therefore be safer in such conditions. Notice 
also that, in terms of coherence preservation and stability, 
the H-protocol shows (as intuitively expected) an inter- 
mediate performance between the A- and R-protocols. 

To summarize: In the logical-IP frame, the effects 
of the RTN can be reduced not only under determinis- 
tic pulses, but also with a randomized control, though 
a comparatively shorter pulse separation is needed in 



A. Free solution for time-independent interaction 
Hamiltonian 

As in Sect. IV, we first focus on understanding the 
controlled dynamics in a frame that explicitly removes 
both the control field and the free evolution due to Hs (S) 
1 + 1 (g) He ■ Let us recall some known results related to 
the uncontrolled dynamics. We have [9I l40j| 

U\t,to) = exp 1^ ® Y.^ble^'^'^'^^kit - to) - h.c.]| , 

(48) 

where 

Cfe(At) = ^(l-e''^^^'). (49) 

Under the standard assumptions that the qubit and the 
environment are initially uncorrelated, 

Ptot(io) = pih) (8) psito) , 
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and that the environment is in thermal equihbrium at 
temperatm'e T [the BoUzmann constant is set =1] 

fc k 

the trace over the environment degrees of freedom may 
be performed analytically, leading to the following ex- 
pression for the qubit coherence: 

pi^{t) - Poi(to)nTrfc{P£,fc(T)I?[e*'^'=*»a(i - ^o)]} 
fe 

= Poi(io)exp[-r(t,to)] . (50) 

Here, I?(^fc) — exp(fej,^fc — bk£X) is the harmonic displace- 
ment operator of the fcth bath mode, and the decoherence 
fimction r(i,to) is explicitly given by 

r(Mo) = E^^^^y^-th(^) . (51) 

k 

In the continuum limit, substituting '^^5{uj — ujk)\gk\'^ 
by the spectral density I{lo), one finds 

r(t, io) - 4 y dui{uj) coth (— j y — m . 

(52) 

For frequencies less than an ultraviolet cutoff Wj,, 
may be assumed to have a power-law behavior, 

/(c.) = . (53) 

The parameter a > quantifies the overall system-bath 
interaction strength and s classifies different environment 
behaviors: s = 1 corresponds to the Ohmic case, s > 1 to 
the super-Ohmic and < s < 1 to the sub-Ohmic case. 

B. Randomly controlled decoherence dynamics: 
Analytical solution and error bound 

Remarkably, the dynamics remains exactly solvable in 
the presence of randomized BB kicks. We focus first on 



the R-protocol viewed in the logical-IP frame. Between 
pulses the evolution is characterized by Eq. H48|) . while 
at tj Eq. applies. Using Eq. (^UJ), the propagator 
in the logical-IP frame, apart from an irrelevant overall 
phase factor, may be finally written as 



U\tMM) = expj^ ®^[6te^"'='"7?f (M, At) -h.c.]| , 

(54) 

where 



M-l 

%«(M, M) = Xje'^^'^'^kiAt) . (55) 

Under the uncorrelated initial conditions specified above 
and thermal equilibrium conditions, the qubit reduced 
density matrix is exactly computed as 



piiitM) = poi{to)llTrk{pEATMe''^'''"VkiM,At)]} 

k 

= Poi(to)e-^"-^*"'*°^ (56) 



Because Xj in Eq. H55|) can be ±1 at random, each ele- 
ment in the sum corresponds to a vector in the complex 
plane with a different orientation at every step At. Thus, 
the displacement operator above may be suggestively in- 
terpreted as a random walk in the complex plane. 

The decoherence function ri^(tM; to) is now given by 



fc ^ ^ 



which, in the continuum limit, becomes 



rii(iAf,to) =4 



du;I{u;) coth (—J ^ '- 



M-l M-j-1 

M + 2J2 cos(j ujAt) J2 ^' Xi+i 



=0 



(58) 



The decoherence behavior under the A-protocol is ob- suit of dcterministically controlled decoherence , which 
tained by letting Xj = (^l)"'- We then recover the re- may be further simplified as [IJ, [43l |. 
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rD(iM,to) =4 



c;w/(u;)coth(^) 



COs[uj{tM - to)] 2 

r LRU 



LoAt 



(59) 



Before proceeding with a numerical comparison be- and the fact that 
tween Eqs. ()58|l - (|59|l . some insight may be gained 

from an analytical lower bound for the average E(x/x/+j) — 1^ "^'^'^^) = j 

E(exp[— rfl(tM, ^o)])- According to Jensen's inequality, 

E(/(x)) > f{E{x)) for any convex function /. Using this we have the following lower bound 



E(exp[-rfl(iM,to)]) >exp[-E(rfl(tM,to))] , 
E{Tn{tM, t,)) = AM r du^Iiu.) coth ( ^) ^ ' '""f^'^ 



(60) 



In Fig. 1^ we compare the coherence decay correspond- 
ing to the absence of control (|52(l . to the A-protocol H59() 
and to the lower bound Ht)U|) . Two limiting cases of high 
and low temperature, T S> Wc and T <C are con- 
sidered. The high temperature limit corresponds to an 
effectively classical bath, where the properties of the en- 
vironment are dominated by thermal fluctuations. In the 
absence of control, decoherence is very fast on the time 
scale determined by the bath correlation time — uj~^, 
hence coherence preservation requires very short intervals 
between pulses. The A-protocol shows the best perfor- 
mance. The actual randomized performance may, how- 
ever, be significantly better than the lower bound in this 
temperature regime, though they never surpass the de- 
terministic case [see next subsection]. 

In the case of low temperature, or fully quantum bath, 
decoherence is much slower and a richer interplay be- 
tween thermal and vacuum fluctuations occurs. Larger 
values of At may then be analyzed before total coher- 
ence loss takes place. The interesting phenomenon of de- 
coherence acceleration 0, Is^ , which may happen when 
ujcAt > 1, may now be observed. For short At, the 
A-protocol is again more efficient, though not signifi- 
cantly better than the lower bound. For large At, pulses 
induce destructive interference and the A-protocol per- 
forms even worse than the lower bound. In such situation 
the best option is simply not to act on the system. 

When LOcAt < 1, some general insight may be gained 
by comparing appropriate limits of the lower bound and 
the deterministic decoherence function. First, by Taylor- 
expanding up to second-order in At we have 

£{rn{tM,to)) ~ (61) 
2(tM - to) At [ du; I{lo) coth (^) , 



whereas 

TD{tMM~ (62) 
At^ dwl{ijj) coth (^) {1 - COs[ijj{tM - to)]} ■ 

Therefore, in the limit of very short At, the lower bound 
approaches the ideal situation of total suppression of de- 
coherence linearly in At, while for the A-protocol this 
occurs quadratically. 

This analysis may be further extended by studying the 
two limits of Eq. I|61|l with respect to temperature. Con- 
sidering the spectral density of Eq. ifCT) , we have 

T » We : £{TR{tMM)) = O {aTul{tM - to)At) , 
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FIG. 9: (color online) Decoherence rate from a bosonic Ohmic 
bath. Here and in the following figures, time is measured in 
units of r"\ a = 0.25, Wc = 100. Left panel: T = IO^lJc and 
cJcAt — 0.1. Right panels: T — 10~^ajc, top: UcAt — 0.1 and 
bottom: cJcAt = 2.5. (Green) Solid line: No control; (blue) 
stars: A-protocol; (red) dashed line: Lower bound, Eq. 1)60^. 
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E(r«(tM,to)) - O {auj'+\tM - to)At) 



Thus, a sufficient condition under which random control 
avoids decoherence is 



aTuj'JtM -to)At <t: 1, 
{tM -to)At<^l 



aoj' 



s+l 



This should be compared with the general bound given 
in Theorem 2 of M. We will see in Sect. V.D that, in 
the physical frame, the qubit frequency too also plays an 
important role. 

A similar analysis for the A-protocol may be effected 
using Eq. For a;c(iM — ^o) ^ 1, the decoherence 

function decays quadratically in time, giving 



T » : ToitM, t„) = O {aTuj^+^tM - At") 



T^Lu,: ToitM, to) = O iaul+^{tM - tofAt") 



while for uJc{tM — to) 3> 1 the decoherence function be- 
comes independent of tM — to, and we get 

T > Lj, : ToitM, to) = O [aTiulAe) , 
T < : VoitM, to) = O [aojl+^At^) . 

These should be compared with the error bound of The- 
orem 3 in [3^. Based on these considerations, random 
pulses may hardly be expected to outperform determinis- 
tic controls in the limits discussed above. Still, it remains 
interesting to quantitatively see what the actual perfor- 
mance is for intermediate At and/or hybrid schemes, for 
instance with respect to acceleration. Moreover, further 
changes may be expected when some time dependence 
exists in the system parameters, for instance in the cou- 
pling strength to the environment [see Sect. V.E]. 



C. Randomly controlled decoherence dynamics: 
Numerical results 

Based on the exact result of Eq. H58|l . we now present 
a comparison of the average decoherence suppression 
achievable by the protocols described in Sect. II. A fixed 
time tf divided into an increasing number M of intervals 
At is considered. 

Fig. EH compares the average E(exp[— ri^(iM, ^o)]) in 
the limit of high temperature, T — lO^cJc, for a system 
evolving under the A- and R-protocols. For the fixed 
times chosen, ujctf = 0.5 (upper panel) and Uctf = 1 
(lower panel), the coherence element has already practi- 
cally disappeared and cannot be seen in the figure, while 
the A-protocol is able to recover it even for very few 
cycles. The values of exp[— rij(tM, ^o)] for different re- 
alizations are widely spread between the worst case cor- 
responding to all A's=0 and the efficient realizations in- 
volving several spin flips. As a consequence, the average 
converges to 1 slowly and has a large standard deviation 
[5^ . Notice, however, that it is significantly better than 
the lower bound. 
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FIG. 10: (color online) Decoherence rate for a high tempera- 
ture Ohmic bath, T = lO^Wc. Upper panel: uJctf = 0.5; lower 
panel: uJctf = 1. (Blue) stars: A-protocol; (black) circles: 
Average over 10^ realizations and respective standard devi- 
ations for the R- protocol [59l| : (black) squares: Expectation 
value (taken over all 2*^ realizations) for the R-protocol; (red) 
dashed line: Lower bound. 



The results in the low-temperature limit, T = 10~^u!c, 
are shown in Fig. ^] Here, thanks to the fact that de- 
coherence is overall slower, longer evolution times may 
be chosen: ujctf = 1 (upper panel) and coctf = 10 (lower 
panel). For the latter choice, in particular, when M < 10, 
decoherence enhancement occurs and, interestingly, the 
results for the A-protocol are worse than those for the 
R-protocol. However, it takes a much smaller At for the 
R-pulses to finally cross the line that separates enhance- 
ment from decoherence reduction. Notice also that the 
values of exp[— rfl(tjv/, ^o)] for different realizations are 
not so spread and the standard deviations are narrower 
than in the high temperature limit. In addition, the aver- 
age over realizations is very close to the lower bound, to 
the point that they cannot be distinguished in the figure. 

We now extend our comparison to the three remain- 
ing protocols of Fig. ^ see Fig. ^1 We choose a high 
temperature bath with coct / = 1 (upper panel) [low tem- 
perature, in this case, leads to similar results], and a low 
temperature bath with LUctf = 10 (lower panel). The S- 
protocol shows the best performance, which is evident in 
the upper panel, but hardly perceptible in the lower one. 
Due to the different rearrangement of the time interval 
between pulses for this protocol, it does not correspond 
to any of the 2^^ realizations of random pulses as consid- 
ered here and represent a special scheme separated from 
the others. The performance of the LS-protocol, which 
has half the number of 7r-pulses used in the A-protocol, 
turns out to be better in all cases of a high temperature 
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FIG. 11: (color online) Decoherence rate for a low tempera- 
ture Ohmic bath, T = lO^^ujc- Upper panel: ujctf = 1; lower 
panel: ujctf = 10. (Green) Solid line: No control; (blue) 
stars: A-protocol; (black) circles: Average over 10^ realiza- 
tions and respective standard deviations for the R-protocol; 
(black) squares: Expectation value for the R-protocol. 



FIG. 12: (color online) Decoherence rate. Upper panel: High 
temperature Ohmic bath, T = lO^tJc, and ujctf — 1. Lower 
panel: Low temperature Ohmic bath, T — 10~^ujc, and fixed 
time uJctf = 10. (Green) Solid line: No control; (blue) stars: 
A-protocol; (purple) plus: H-protocol; (orange) diamonds: 
LS-scheme; and (black) dot-dashed line: S-protocol. Aver- 
age for the H-protocol is taken over 10'^ realizations. 



bath, but worse for a fully quantum bath with large tOctf. 
This explains why the H-protocol, which combines sym- 
metrization and randomness, also performs better than 
the A- in the high temperature limit. 

To summarize: In terms of performance, we have in 
decreasing order: S-, LS-, H-, A-, R-protocols for high 
temperature; and S-, A-, H-, LS-, R-protocols for low 
temperature once the number of pulses are sufficient to 
start slowing down decoherence. Different protocols be- 
come again, as expected, essentially equivalent in the 
limit A< — > 0. For finite pulse separations, in the consid- 
ered case of a time-independent Hamiltonian, it is always 
possible to identify a deterministic protocol showing the 
best performance. However, if a balance is sought be- 
tween good performance and protocols minimizing the 
required number of pulses, then the H-protocol again 
emerges as an interesting compromise. Note, in particu- 
lar, that the latter outperforms the standard A-protocol 
in some parameter regimes. 



D. Randomized decoupling in the physical frame 

We now investigate under which conditions decoher- 
ence suppression is attainable in the physical frame, when 
the system is subjected to randomized control. Because, 
in this frame, the qubit natural frequency plays an impor- 
tant role, random decoupling also depends on how small 
At can be made with respect to tq = WjT^. 



The reduced density matrix is obtained following the 
same steps described so far, but in order to retain the ef- 
fects of the system Hamiltonian, the transformation into 
the interaction picture is now done with respect to the 
environment Hamiltonian only - hence the superscript 
IE. Upon tracing over the environment degrees of free- 
dom, we are left with the reduced density operator in the 
Schrodinger picture. 

The unitary operator between pulses is 



U'^{tj+i,t,) - exp (^^i'^a.At 



(63) 



expjy®^[&(At)6te-^*^-h., 



while at tj it is given by 



P. 



IE 



exp 



(64) 



By additionally moving to the logical frame we get 

M-l 



exp 



^[6le-'=*«r;f(M,Ai) 



(65) 



Tracing over the environment and taking the expectation 
over control realizations leads to the coherence element 
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in the logical frame: 



E 

fe=i 



where .n,'''^'' 



(66) 

Af-i is given by Eq. (gj). Thus, in addition to the 
decoherence described as before by Eq. H58|) , we now have 
ensemble dephasing due to the fact that each realization 
carries a different phase factor proportional to loq- 
The results for the ratios 



Fi(tAf) 



E(/5oi(^A/)) 
Poi^o) 



and ¥2(1 m) 



E(/5oifa/)) 
Poi(io) 



(67) 

in the logical and logical-IP frames for the system, re- 
spectively, are summarized in Fig. 1131 where At is fixed 
and the system is observed at different times. Both a high 
temperature and a low temperature scenario are consid- 
ered. The phase for each realization in the logical frame 
is mostly irrelevant when ujq <ti 1/A.t. The outcomes 
of the average over all realizations in both frames are 
then comparable, independently of the bath temperature. 
The situation changes dramatically when the spin-flip en- 
ergy becomes large, the worst scenario corresponding to 
kn/{2At), with k odd. Here, because S^'^'' 



Wo 



even number, each realization makes a positive or a neg- 
ative contribution to the average, which may therefore 
be very much reduced. Such destructive quantum inter- 
ference is strongly dependent on the bath temperature. 

Among all random pulse realizations, the most effec- 
tive at suppressing decoherence are those belonging to 
the smaller ensemble of the H-protocol. None of them 
carries a phase, so they always make large positive con- 
tributions to the total ensemble average. In a high tem- 
perature bath, the realizations that can make negative 
contributions have often tiny values of e~'"^*^*"'*°-', which 
explains why even in the extreme case of ujq = fc7r/(2At) 
the R-protocol can still lead to some decoherence reduc- 
tion. In a low temperature bath, on the other hand, 
decoherence is slower and for the time considered here, 
the values of e~^"*^**^'*°' for all realizations are very close, 
which justifies their cancellation when loq = kn / {2\t). 

In the physical frame, the average for the density ma- 
trix depends on the initial state of the system as 



E(poi(^m)) 



/9oi(io) + Pio(io) 



2 • 2 



M-l 



E cos[si^)c.oA<]e-^«'(*--*«) . 

k=l 

As already discussed in Sect. Ill, the problem associated 
with population inversion may be avoided if a classical 
register is used to record the actual number of spin flips. 

To summarize: Two conditions need to be satisfied for 
the R-protocol to become efficient in reducing decoher- 
ence: WcAt <C 1 and also woAt ^ 1. Notice, however, 
that when randomness and determinism are combined in 



a more elaborated protocol, such as the H-protocol, no 
destructive interference due to ujq occurs. In addition, 
the hybrid scheme is still capable to outperform the A- 
protocol in appropriate regimes. 



E. Time-dependent coupling Hamiltonian 

As a final example, imagine that the coupling param- 
eters gk{t) between the system and the environment are 
time dependent and let us for simplicity work again in 
the logical-IP frame. The total Hamiltonian is given by 
Eqs. Q-© with 0^0 (t) = ijJQ and ^ = 1. Two illustrative 
situations are considered: gk{t) changes sign after certain 
time intervals, or it periodically oscillates in time. 



1. Instantaneous sign changes 
Suppose that gk{t) = gkD{t), where 
^(t) = (-l)L^J 



(68) 



describes instantaneous sign changes of the coupling pa- 
rameter after every interval 3/(10a;c)- For a high tem- 
perature bath and a fixed time tf — 1/lUc, Fig. 1141 shows 




FIG. 13: (color online) Ratios Fi{t) and F2{t) in the logi- 
cal and logical-IP frames, respectively. A fixed time inter- 
val At = 1/(10lJc) is taken. Upper panel: High tempera- 
ture Ohmic bath, T — W^loc, lower panel: Low temperature 
Ohmic bath, T — 10^^ ujc- (Green) Solid line: No control; 
(blue) stars: A-protocol; (purple) plus: H-protocol; (black) 
squares: R-protocol in the logical-IP frame; (red) up trian- 
gles: R-protocol in the logical frame with small frequency 
LOoAt = 10~^; (red) down triangles: R-protocol in the logical 
frame with large frequency ujoAt — n/2. Average performed 
over all realizations. 
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FIG. 14: (color online) Decoherence rate for a high tempera- 
ture Ohmic reservoir, T = W^lOc, with alternating couplings. 
(Blue) stars: A-protocol; (black) circles: R-protocol; (pur- 
ple) plus: H- protocol. The interval considered is tf = 1/ujc- 
Averages taken over all possible realizations. 
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that the results for the A-protocol exhibit a drastic drop 
when At = tf/4: and At = tf/10. This is due to the fact 
that some of sign changes happen very close to or coin- 
cide with some of 7r-pulses of the deterministic sequence, 
canceling their effect. In contrast, the occurrence of spin 
flips in randomized schemes is irregular, so that the lat- 
ter are more protected against such "resonances" and 
steadily recover coherence as At decreases, even though 
at a slower pace. Note that when dealing with the S- or 
LS-protocols, the same sort of recoil should be expected 
for different time dependences and different values of At. 



2. Periodic modulation 

Assume that the coupling parameter is given by 
9k{t) = 9kG{t), where G{t) = cos{p nujct) sm{q nujct) and 
\p — q\ is small. This function has two superposed pe- 
riodic behaviors, one with a long period and the other 
fast oscillating. The fast oscillations are shown in the 
left upper panel of Fig. ^1 

We consider a high temperature bath, T = lO^Wc- The 
right upper panel shows the qubit decoherence in the ab- 
sence of pulses. The oscillations in the decay rate are 
related to the oscillations in the interaction strength be- 
tween the system and the bath. In the lower panel, we fix 
a time tf = l/ujc and compare the decoherence rate for 
the cases of absence of control, A-, H- and R-protocols. 
When At = tf/Q the result for the A-protocol suddenly 
becomes even worse than not acting on the system. Ran- 
dom pulses, on the contrary, do not show any significant 
recoil. The reason for the inefficiency of the A-protocol 
when M — 6 becomes evident from the left upper panel 
of Fig. El Vertical dashed lines indicate where the pulses 
occur. They mostly coincide with the instants where G{t) 
also changes sign. For the LS-protocol, similar unfavor- 
able circumstances happen for different values of At and 
similar behaviors should be expected for other determin- 
istic protocols and functions gk{t). 



FIG. 15: (color online) Upper left panel: Function G{t) = 
cos{pTTuict)sin{qnujct), for p = 2.95 and q — 3.25. Upper 
right panel: decoherence rate in the absence of control. Lower 
panel: Decoherence rate for a fixed time interval LOctf = 1. 
A high temperature Ohmic bath, T = lO'^Uc, is considered. 
(Green) Solid line: No control; (blue) stars: A-protocol; 
(black) circles: R-protocol; (purple) plus: H-protocol. Av- 
erages taken over 10'^ realizations. Standard deviations for 
the R-protocol are shown. 

To summarize: The above examples again reinforce 
the idea of enhanced stability of randomized controls, 
and suggest that randomization might represent a safer 
alternative in reducing decoherence when limited knowl- 
edge about the system-bath interaction is available. 



VI. CONCLUSION 
A. Summary 

A quantitative comparison between deterministic and 
randomized control for the most elementary target sys- 
tem, consisting of a single (isolated or open) qubit, was 
developed in different frames. The main conclusions 
emerging from this study may be summarized as follows. 

First, it is always possible to identify conditions un- 
der which purely random or hybrid schemes succeed at 
achieving the desired level of dynamical control. Frame 
considerations play an important role in specifying such 
conditions, satisfactory performance in a given frame be- 
ing ultimately determined by a hierarchy of time scales 
associated with all the dynamical components in the rel- 
evant Hamiltonian. While all protocols become essen- 
tially equivalent in the limit of arbitrarily fast control, 
the behavior for finite pulse separation is rich and rather 
sensitive to the details of the underlying dynamics. As a 
drawback of pure random design, ensemble average tends 
to introduce, in general, additional phase damping, which 
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may be however circumvented by combining determinism 
and randomness within a hybrid design. 

Second, for time-independent control settings in this 
simple system, it was always possible to identify a deter- 
ministic protocol with best performance. While deter- 
ministic schemes ensuring accurate averaging of a known 
interaction always exist in principle [l|, such a conclu- 
sion remains to be verified under more general circum- 
stances, in particular access to a restricted set of control 
operations. The hybrid protocol proved superior to the 
pure random schemes, as well as to standard asymmetric 
schemes in certain situations. 

Third, for time- varying systems, randomized protocols 
typically allow for enhanced stability against parameter 
variations, which may severely hamper the performance 
of deterministic schemes. Pure random design tends to 
perform better than hybrid in this respect, both choices, 
however, improving over purely cyclic controls under ap- 
propriate conditions. 

Overall, hybrid design emerges as a preferred strategy 
for merging advantageous features from different proto- 
cols, thereby allowing to better compromise between con- 
flicting needs. 



B. Outlook 

From a conceptual standpoint, it is intriguing to realize 
that complete suppression of decoherence remains pos- 
sible, in principle, by purposefully introducing a proba- 
bilistic component in the underlying control, and perhaps 
surprising to identify cases where this leads to improved 
efficiency over pure deterministic methods. 

In a broader context, however, it is worth mentioning 
that the philosophy of recognizing a beneficial role of ran- 
domness in physical processes has a long history. Within 
NMR, the stochastic averaging of intermolecular interac- 
tions in gases and isotropic liquids due to random transla- 
tional and re-orientational motions may be thought of as 
a naturally occurring random self-decoupling process . 
In spectroscopic applications of so-called stochastic NMR 
and stochastic magnetic- resonance imaging [60| , spin ex- 
citation via trains of weak RF pulses randomly modu- 
lated in amplitude, phase, and/or frequency are used to 
enhance decoupling efRciencies over a broader frequency 
bandwidth than attainable otherwise. Even more gen- 
erally, the phenomenon of stochastic resonance |6ll l62j| 
is paradigmatic in terms of pointing to a construc- 
tive role of noise in the transmission of physical sig- 
nals. Within QIP, strategies aimed at taking advan- 
tage of noise and/or stochasticity have been considered 
in contexts ranging from quantum games '63l, to quan- 
tum walks .64], dissipation-assisted quantum computa- 
tion 65], as well as specific coherent-control |^ and 
quantum simulation |67l] scenarios. Yet another sug- 
gestive example is offered by the work of Prosen and 
Znidaric, who have shown how static perturbations char- 
acterizing faulty gates may enhance the stability of quan- 



tum algorithms 68]. More recently, as mentioned, both 
pure random ^3^] and hybrid [s^ active compensation 
schemes for static coherent errors have been proposed. 
While it is important to stress that none of the above 
applications stem from a general control-theoretic frame- 
work as developed in |36j , it is still rewarding to fit such 
different examples within a unifying perspective. 

Our present analysis should be regarded as a first step 
toward a better understanding and exploitation of the 
possibilities afforded by randomization for coherent and 
decoherent error control. As such, it should be expanded 
in several directions, including more realistic control sys- 
tems and settings, and fault-tolerance considerations. 
While we plan to report on that elsewhere, it is our hope 
that our work will stimulate fresh perspectives on fur- 
ther probing the interplay between the field of coherent 
quantum control and the world of randomness. 
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APPENDIX A: CONTROL HAMILTONIAN 

The control Hamiltonian is designed according to the 
intended modification of the target dynamics in a desired 
frame. Throughout this work, our goal has been to freeze 
the system evolution by removing any phase accumulated 
due to the unitary evolution, as well as avoiding nonuni- 
tary ensemble dephasing and decoherence. As clarified 
below, this requires the use of identical n-pulses in the 
physical frame. This condition may be relaxed at the 
expense of no longer refocusing the unitary evolution. 

Let the control of the system be achieved via the ap- 
plication of an external alternating field (e.g., a radiofre- 
quency magnetic field), 

H,{t) = J2 (t) = ^''^ ^o^t'^* + (^)]'^- ' 

j 3 

(Al) 

where the carrier is tuned on resonance with the qubit 
central frequency, to = ujq. As described in the text, 
yO) (t) = y [e(t - tj) -Q{t-tj- t)] and each pulse hap- 
pens at tj, having duration t and amplitude V. Upon 
invoking the rotating wave approximation (RWA) , hence 
neglecting the counter-rotating terms (T_ exp[— i(woi + 
(/7j(i))] and (t+ exp[i(wot + V'j (*))] i^i I|A1() . the control 
Hamiltonian given in the main text is found. The func- 
tion ifj (t) characterizes the phase properties of the pulses 
we deal with. We compare two relevant possibilities: 
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(i) (pj{t) = —ujQtj for each j: This means that the pulses 
are identical in the physical frame, as used in this work. 

(ii) (pj{t) — for all j: This means that the pulses are 
identical in the physical frame only if separated in time 
by a multiple of 27r. 

The propagator corresponding to the above choices 
may be in general obtained by seeking a transformation 
which removes the time dependence of Hc^\t) within 
each pulse. A transformation to an absolute frame ro- 
tating with the carrier frequency, which on resonance is 
identical with the interaction picture, leads to 



ijfW(i) = exp[iwotCTj2]ff^^')(t)exp[-iwotcrz/2 
= V^O)(t)expU^Ja,exp 



Thus, the choice ^j{t) — corresponds to pulses which 
are translationally invariant in time in this frame. In case 
(i), the above transformation does not remove time de- 
pendence, which would instead be accomplished by mov- 
ing to a relative rotating frame via a rotation Uz(t — tj) = 
exp[iujo{t — tj)az/'2]. From the above expression, the 
interaction-picture propagators for an instantaneous tt- 
pulse applied at tj are found, respectively, as 



(i) P/ = exp 

(ii) Pj — exp 



—i—a., 
2 



exp 



TT 

-i—a^ 
2 



exp 



.Wot, 



relation 

Pj = exp[—iuJotj (7 z /'2]Pj exp[ia;oijO'z/2] 
leading to the propagators 



(i) Pj = exp 

(ii) Pj ~ exp 



TT 

. t^otj 



(A2) 









. ^otj 


exp 


. 2 


exp 


[ 2 "^J 



Thus, pulses with ipj{t) = —LOotj are confirmed to be 
translationally invariant in time in the physical frame, as 
directly clear from the dependence uJo{t — tj) in HA1|) . 

The difference between the two choices to the con- 
trol purposes becomes evident by considering the A- 
protocol on the isolated qubit. From HA2|I . the propaga- 
tors U{t2,to) = P2Uo{t2,ti)PiUo{ti,to) in the physical 
frame are, respectively. 



1, 



(i) U{t2,to) 
(u) U{t2, to) = - exp [-iujoit2 ~ ti)az] , 



(A3) 



We can then return to the Schrodinger picture using the 



which leads to the conclusion that refocusing in the physi- 
cal frame may only be achieved with identical pulses, that 
is, if (pj{t) = ujotj. Clearly, for the choice (pj{t) = 0, the 
accumulated phase may only be disregarded in the frame 
rotating with the frequency uq. Both choices are equally 
useful if decoherence suppression becomes the primary 
objective in the open system case. 
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